The meandering instability of a viscous thread 
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A viscous thread falling from a nozzle onto a surface exhibits the famous rope-coiling effect, in 
which the thread buckles to form loops. If the surface is replaced by a belt moving with speed U, 
the rotational symmetry of the buckling instability is broken and a wealth of interesting states are 
observed [See S. Chiu- Webster and J. R. Lister, J. Fluid Mech., 569, 89 (2006)]. We experimentally 
studied this "fluid mechanical sewing machine" in a new, more precise apparatus. As U is reduced, 
the steady catenary thread bifurcates into a meandering state in which the thread displacements 
are only transverse to the motion of the belt. We measured the amplitude and frequency uj of 
the meandering close to the bifurcation. For smaller U, single-frequency meandering bifurcates to a 
two-frequency "figure eight" state, which contains a significant 2iL} component and parallel as well as 
transverse displacements. This eventually reverts to single-frequency coiling at still smaller U. More 
complex, highly hysteretic states with additional frequencies are observed for larger nozzle heights. 
We propose to understand this zoology in terms of the generic amplitude equations appropriate for 
resonant interactions between two oscillatory modes with frequencies cj and 2uj. The form of the 
amplitude equations captures both the axisymmetry of the U = coiling state and the symmetry- 
breaking effects induced by the moving belt. 

PACS numbers: 82.40.Bj, 47.20. Gv, 47.20. Ky 



The fascinating instabilities of very viscous fluids are 
familiar to anyone who has poured syrup onto a pan- 
cake [1 . The thread of syrup elongates as it falls, and 
exhibits a buckling instability as nears the surface [5]. 
The buckling is the result of a competition between bend- 
ing and axial compression, and causes the thread to loop 
itself into coils. This "liquid rope coiling effect" [3] has a 
long history comprising almost a half century of experi- 
ments m H S 13 El [7] , scaling analyses [3 H IH] , analytic 
theory and numerical simulation [7] [TOl HU [12 [SI Ej . 
It has been observed that the coils tend to pile up on 
the surface, forming a tall column [31 HI [S] and that the 
frequency of coiling exhibits a complex, multivalued de- 
pendence on the fall height [5 , 7\. Here we consider a vari- 
ation on this classic experiment, the case when the liquid 
thread falls onto a moving surface, a situation known as 
the "fluid-mechanical sewing machine" [TS] owing to the 
variety of patterns deposited on the surface. 

The moving surface breaks one of the basic symmetries 
of the "pure" rope coiling problem, and has the effect of 
unfolding the coiling instability into a rich panoply of 
distinct bifurcations. This arrangement is also very con- 
venient experimentally; the moving surface is provided 
by a belt which is cleared of fluid before it returns, elimi- 
nating the pile- up of viscous fluid at the point of contact. 

The main experimental control parameters are the 
height H of the nozzle from which the thread descends, 
and the speed of the moving surface U. Less important 
parameters, which we hold fixed in this work, are the 
volumetric flow rate of the liquid Q and the diameter of 
the nozzle d. The Newtonian fluid is characterized by its 
density p, kinematic viscosity u and surface tension a. 
Non- Newtonian effects are negligible. 



We present an experimental study of the "stitch" pat- 
terns made by the thread as it is laid down on the belt, 
as a function of H and U. We flrst discuss a survey of 
the states in the H —U plane, with denser coverage than 
was possible in Ref. [TS]. Then we focus on the simplest 
bifurcation, that from the straight catenary state to me- 
andering. We show that this state is well described as 
a forward Hopf bifurcation and compare the onset belt 
speed and Hopf frequency to recent linear stability the- 
ory [16'. We examine the nonlinear saturation of the 
meandering amplitude and qualitatively explain it with 
a simple kinematic model of an inextensible thread. We 
then essay the task of understanding the more complex 
"figure 8" and other patterns that appear as the belt 
speed is lowered. We propose a general amplitude equa- 
tion approach, based on the idea that the motion of the 
belt can be treated as a perturbation to the pure rope 
coiling problem which has 0(2) symmetry. We propose 
a set of coupled amplitude equations which minimally 
break 0(2) symmetry, and mix amplitudes for rotating 
wave modes with frequencies uj and 2a;. We conjecture 
that all the complicated stitch patterns can be captured 
by this approach, which could be generalized to include 
more modes as they become excited near the multifre- 
quency regime. These modes can be directly measured 
by visualizing longitudinal and transverse motions of the 
thread by viewing from the side. 

This paper is organized in the following way. In section 
|T] we discuss the fluid properties and the experimental 
apparatus. In section ITT] we make comparisons between 
our parameters and the corresponding zero belt speed, 
pure coiling problem, before turning to the state diagram 
in section jll A| and the meandering threshold in sections 



II B and II C In section [TlD| we discuss experimental ob- 
servations of the states beyond meandering. Section III 
presents our general amplitude formalism which provides 
a consistent framework for interpreting the experimen- 
tal results. We have placed a detailed discussion of the 
derivation of the amplitude equations in an Appendix. 
Section IIVI is a brief conclusion. 



I. EXPERIMENT 

Our experimental apparatus was a redesigned version 
of the one described in Ref. [TF . It is shown schematically 
in Fig. IT] The main differences were improved control 
and measurement of the belt speed, and the use of very 
stable silicone oil as the working fluid. 

Silicone oil offers several advantage over the Lyle's 
Golden Syrup used in previous experiments |15j . Al- 
though its viscosity cannot be varied by adding water 
as described in Ref. [15], silicone oil is not susceptible to 
evaporation and its viscosity varies much less with tem- 
perature than that of sugar syrup. This makes it possible 
to do highly reproducible experiments over a long period 
of time. It also has about one third the surface tension of 
syrup, which reduces the thread-thread and thread-belt 
interaction after the fluid is laid down on the belt. 

We used Dow Corning 200© fluid, which has a nominal 
kinematic viscosity of 30,000 cSt. To fully characterize 
the fluid, we measured its viscosity and the temperature 
coefficient of its viscosity using a TA ARIOOO rheome- 
ter. The same instrument was used to confirm that the 
fluid is Newtonian to an excellent approximation. We 
also measured the fluid density and its temperature co- 
efficient using an Anton-Parr DMA 5000 densitometer. 
The results of these measurements are shown in Table H] 

The temperature of the experiment was recorded dur- 
ing each run, but not otherwise controlled. The only 
significant variation in the fiuid properties was therefore 
due to the temperature, and not to the degradation of 
the fluid. Averaged over all experiments, the mean tem- 
perature was 22.2±0.8°C. Accounting for this systematic 
error, the fluid had a density of p = 996 ± 8 kg/m'^ and 
a kinematic viscosity of ly — (2.77 ± 0.02) x lO^'^ m^/s, 
or 27, 700 ±200 cSt. 

We used an internally reinforced, toothed timing belt 
of width 16 mm. The belt was looped over toothed pul- 
leys and driven with a speed-controlled DC motor. The 
speed of the belt U was measured using an averaging stor- 
age oscilloscope to record the intensity of a laser beam 
that was interrupted by the teeth of the belt. The fre- 
quency of the passage of the teeth, measured by an in- 
ternal function of the oscilloscope, could be used to find 
U to within a few % . The reinforced timing belt resisted 
stretching or slipping, even under considerable pressure 
from the plastic scraper used to remove the silicone oil. 
After scraping, the belt had only a very thin wetting layer 
of oil left on it. The excess oil dripped off the scraper into 
a beaker, and could be recirculated. 



Density p 


996 ± 8 kg/m^* 


Temperature coefficient of p 


-0.0885 %/°C 


Molecular viscosity rj 


27.63 ± 0.05 Pa s 


Temperature coefficient of rj 


-1.6 %/°C 


Kinematic viscosity u = rj/ p 


(2.77 ±0.02) X 10"^ mVs 


Maxwell time tm 


< 1.8 ms 


Young's modulus G 


< 14 kPa 


Surface tension a 


21.5 mN/m 



TABLE I: The measured properties of the silicone oil at the 
mean experimental temperature of 22.2 ± 0.8° C. The size of 
the measured Maxwell time and Young's modulus show that 
viscoelastic effects are negligible. The surface tension is taken 
from the manufacturer's data sheet. 



In all the experiments reported here, we used a noz- 
zle with d = 8.00 ± 0.02 mm to deliver a fixed, constant 
volumetric fiux of oil Q. The oil was supplied by a reser- 
voir and driven by a gravity feed system with a fixed 
42 cm head. The reservoir and the nozzle were attached 
to a frame mounted on a screw so that they travelled to- 
gether, maintaining a constant head as the fall height H 
was varied. H was measured with a scale to ±0.05 cm. 

Before and after each run, we weighed the quantity of 
oil delivered in a certain time interval, typically 2 min, 
timed with a stopwatch. We found the mass flux pQ = 
0.0270 ±0.0006 g/s remained constant within experimen- 
tal error, even for runs performed many days apart. All 
the experiments we discuss here used this value of pQ. 

An experimental run consisted of fixing the nozzle 
height H and varying the belt speed U, starting from 
high values of U for which the thread is dragged hori- 
zontally by the belt to form a steady, stretched catenary. 
Lowering the belt speed to zero, and then increasing it 
again produces a sequence of bifurcations, and a sequence 
of "stitch patterns" of the thread on the belt. Care was 
taken to change U slowly near the onset of the first insta- 
bility, to meandering, where the state of the thread re- 
sponded very slowly to small changes in U. We recorded 
the oil pattern on the moving belt from above using a 3 
megapixel digital camera. To avoid heating the oil, and 
to freeze the motion of the thread, the only lighting used 
was a fast photo flash. Typical photos of the main states 
of interest are shown in Fig. [2] These and many other 
states are shown and described in Ref. [TS]. 



II. RESULTS AND DISCUSSION 

In this section we will tour the zoology of states and 
bifurcations exhibited by the viscous thread falling on 
the moving belt. After a general overview, we focus on 
the simplest transition first, the onset of the meander- 
ing state from the stretched catenary as the belt speed 
is lowered. We compare our results to the existing linear 
stability theory for this transition [TB]. Moving into the 
nonlinear regime, we find that the amplitude of the me- 




Scraper 



Toothed belt 



FIG. 1: A schematic of the apparatus. A sihcone oil stream 
with volumetric flow rate Q falls a distance H from a nozzle of 
diameter d onto a belt moving at speed U . The oil is removed 
by the scraper. The laser, which was interrupted by the teeth 
of the belt, was used to measure the speed U . 
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FIG. 2: Photos of some of the states of the thread on the belt. 
The nozzle is on the left and the belt moves from left to right. 
The internal reinforcing mesh of the semi-transparent belt is 
visible, but the surface of the belt is smooth. The scale is 
in mm. (a) The catenary steady state, (b) The meandering 
state, (c) The figure eight state, (d) The translated coiling 
state. 



andering is well described as a forward Hopf bifurcation. 
We then show that the saturation coefficient of the 
amplitude can be deduced reasonably well from a sim- 
ple model of an inextensible thread. Finally, we consider 
the subsequent bifurcations for small nozzle heights H, 
namely the appearance of the figure eight and translated 



coiling states. These and the other more intricate states 
for larger H are probably best understood from the dy- 
namics of their amplitudes. We present a brief outline of 
a bifurcation theoretic approach to this problem. 

At the outset, it is instructive to identify the states of 
"pure" rope coiling that we would expect to see if the 
speed of the belt in our experiments were zero. Theory 
and laboratory experiments show [71 [T31 [T3], that pure 
coiling can occur in four different regimes - viscous (V), 
gravitational (G), inertio-gravitational (IG), and inertial 
(I) - depending on the relative magnitudes of the different 
forces acting on the thread. In the V, G, and I regimes, 
the coiling frequency fi is a single- valued function of the 
fall height H . In IG coiling, by contrast, Q.{H) is multi- 
valued, reflecting the fact that the nearly vertical upper 
part of the thread acts as a distributed pendulum with 
multiple resonantly excited eigenmodes. 

In a typical experiment where the fluid properties, the 
flow rate Q and the nozzle diameter d are held constant, 
the regime that obtains is determined by the fall height 
H. As H increases, the regimes always succeed each 
other in the order V-G-IG-I. However, while the G and I 
regimes can be observed for nearly any choice of Qi d, and 
the working fluid, the V and IG regimes are confined to 
more restricted regions of parameter space. For a given 
experiment, this space is completely characterized by the 
three dimensionless groups. 



Hi 

n2 

Ha 



,5 \ 1/5 



9Q' 



puQ 



1/4 



[610 ±10] 


(1) 


[0.370 ± 0.002] 


(2) 


1.84 ±0.06] 


(3) 



where g is the acceleration due to gravity, a is the surface 
tension, and the values of each group for our experiments 
are given in square brackets. For a fluid such as sili- 
cone oil with a relatively low surface tension, the value 
of Ha has only a minor effect on the coiling behavior, 
which is controlled primarily by Hi and 112. It turns out 
that our values Hi « 610 and 112 ~ 0.37 are within the 
domain where IG coiling is expected to occur for some 
range of heights (See Ref. [7], Fig. 3.) By contrast, the V 
regime does not occur, because the characteristic height 
iyQI 9)^ ~ 3 mm over which gravitational stretching of 
the thread becomes important is less than the diameter 
of the nozzle (as usual with a gravity feed system) . A di- 
rect numerical calculation using the method of [13^ shows 
that G coiling is expected in our experiments for 2 cm 
< H < 5.7 cm, IG coiling for 5.7 cm < H < 8.5 cm, and 
I coiling for H > 8.5 cm. The experiments reported here 
explore the range 2 cm < H < 8 cm, with belt speeds up 
to C/ = 10 cm/s. 
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FIG. 3: The state diagram, as a function of the belt speed U and the nozzle height H. Some of the simpler states are shown 
in the previous figure; more complex states are described in Ref. [15] . 



A. The state diagram 

Fig. [3] shows the states observed in the U-H plane. 
The diagram was made by scrutinizing over 1100 indi- 
vidual photos. In general, the diagram agrees very well 
with Fig. 6 of Ref. [K], although certain details differ at 
larger H. At large U, only the stretched catenary steady 
state is seen. As U is reduced at fixed H, the first bifur- 
cation is to the meandering state. This transition, which 
is non-hysteretic for H < 7.0 cm, is discussed in detail 
below. Above H = 7.0 cm, the meandering state shows 
clear multifrequency behavior and a small hysteresis at 
onset. For H < 5.5 cm, as U is reduced, the meandering 
state switches to a figure eight mode, before making a 
second transition to translated coiling. All these tran- 
sitions are very nearly non-hysteretic. It is not clear if 
the figure eight window extends all the way down to the 
onset of any instability, just below H = 2.0 cm. For 
5.7 cm < H < 7.0 cm, the figure eight state disappears, 
and there is a small hysteretic region where meander- 
ing and coiling co-exist. For H > 7.0 cm, many com- 
plex states are observed, with windows of various pat- 



terns separated by disordered, possibly chaotic regions. 
In contrast to Ref. [TS] , almost no region of slanted loops 
is observed. This may be due to the reduced surface ten- 
sion of the silicone oil, which reduces the tendency of the 
thread to stick to itself. Regions of W-states and dou- 
ble coiling are seen, as well as small re-entrant regions of 
the simpler figure eight and single-frequency meandering 
modes. The latter was not observed in Ref. jTS]. We do 
not observe side-kicks, which may require still larger H. 
A complete elucidation of the multifrequency states for 
H > 7.0 cm must await more precise and better instru- 
mented experiments. For the remainder of this paper, we 
will focus on the simpler bifurcations between the steady 
and meandering states, and on the region of the figure 
eight window for H < 5.5 cm. 



B. The threshold of single frequency meandering 

We now consider the bifurcation from the steady 
stretched state to meandering, which has been analyzed 
at the linear stability level in Ref. [16 . In particular, lin- 



ear theory provides clear predictions for the critical belt 
speed Uc at which meandering begins as U is decreased, 
and for the critical frequency cOc of the meandering state 
which emerges. Quite remarkably, the onset angular fre- 
quency ujc (which, as we will show, is the Hopf frequency 
that obtains for a zero amplitude of meandering) is pre- 
dicted to be nearly identical to the frequency expected 
for pure, nonzero amplitude rope coiling onto a station- 
ary surface for the same H [IB] . 

To accurately locate the onset of meandering, we pho- 
tographed the state while very slowly decreasing and 
increasing U through Uc- We then found the peak-to- 
peak amplitude 2\A\ and wavelength A of the meandering 
thread directly from the digital photos. It was typically 
possible to use data spanning 5 to 10 wavelengths, and 
the resolution of the images was sufficiently high that the 
width of the thread was approximately 10 pixels. The 
smallest amplitudes we could measure were a fraction of 
the width of the thread, while the largest spanned the 
belt. 

We fit the amplitude data to a phcnomcnological am- 
plitude model of the Landau form 



TA = eA- ^J,A\A\'^ + h.o.t. 



(4) 



with e — {Uc — U)/Uc- r is a time scale related to the 
linear growth rate. We consider the steady state of this 
equation, so tA = and use the fit to determine Uc 
and the saturation coefficient ji. Fig. |4] shows such a fit. 
The fits were best near the midrange of the heights H , 
declining in quality for small H due to sparsity of data 
and for larger H due to the growth of the higher order 
terms. Near H fa b cm, the amplitude has an essentially 
perfect square-root dependence \A\ = yje/^i, as shown in 
Fig. |4] For H > 7.0 cm, fitting was no longer possible, 
and we simply estimated Uc by bracketting the onset of 
multifrequency meandering. 

Fig. [5] shows the critical belt speed Uc as a. function 
of the nozzle height H, compared to the prediction of 
the linear stability theory [15]. There are no adjustable 
parameters in the theory. The backward looping portions 
of the theoretical curve mark the inertial-gravitational 
regime where multifrequency coiling would exist in the 
pure rope coiling case. Since our experimental protocol 
was to fix H and vary U, we would expect to observe 
the upper envelope of this curve. The agreement is very 
satisfactory. 

We have no immediate explanation for the slight sys- 
tematic discrepancy by which the theory overestimates 
Uc by an amount that increases with H. Presumably it 
is due to some unaccounted for physical effect, such as 
air drag or surface tension effects at the contact point be- 
tween the thread and the belt. For large H, the thread 
is very sensitive to air currents and the rapidly moving 
belt inevitably disturbs the air nearby. 

The meander frequency to can also be easily obtained 
from the measurements of the wavelength A of the mean- 
der pattern on the belt. The angular frequency is simply 
u! = 2TrU/X. We find that lo decreases linearly with U, 
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FIG. 4; (Color online) The amplitude \A\ of the meandering 
motion as a function of the belt speed U, for H — 5.3 cm. The 
upward (downward) pointing triangles indicate data taken 
with increasing (decreasing) U. The solid line is a fit to Eqn.E] 
with Uc ~ 4.01 cm/s and fi — 4.62 cm~^. No higher order 
terms were required. The smallest nonzero amplitudes here 
are about half the width of the thread. 



so that it is straightforward to extrapolate to obtain lOc 
as shown in Fig. [6j The functional forms of the ampli- 
tude and frequency of meandering shown in Figs. |4] and [6] 
demonstrate that the onset of single frequency meander- 
ing is very well described as a forward Hopf bifurcation 
with Hopf frequency tOc- This onset frequency is also 
predicted by linear stability theory |Tnj; Fig. [t] shows the 
Hopf frequency as a function of H. Here the agreement 
with the theory is excellent, and again there are no ad- 
justable parameters. 

It is interesting that the frequency measurements track 
the lower branch of the multifrequency regime, right to 
the end of its existence. This is a consequence of the ex- 
perimental protocol in which H is fixed. Presumably the 
other frequency branches could be accessed by varying 
both H and U, or perhaps by making finite perturbations 
to the state of the thread. The frequency branches have 
been observed in pure rope coiling experiments O [7j. 
The lower branch and the observation of single frequency 
meandering both end near H — 7.0 cm. Above this value, 
we find mostly multifrequency meandering at onset, as 
shown in Fig. [3] 



C. Meandering amplitude saturation: 
a simple model 

The fits to the amplitude of meandering also give a 
quantity that is not predicted by linear theory, namely 
the coefficient /i of the cubic nonlinear term in Eqn. [4] 
Fig.[8]shows the experimental result. We find that /i is a 
decreasing function of H. The coefficient /x controls the 
dependence of the saturated amplitude on the belt speed. 
It can be predicted by the following argument, which is 
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FIG. 5: (Color online) The critical belt speed Uc vs H, ex- 
tracted from fits to amplitudes like the one shown in the previ- 
ous figure. The solid line is calculated from the linear stability 
analysis of Ref. [TS], with no adjustable parameters. 
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FIG. 6: (Color online) The angular frequency oj of the mean- 
dering motion (open circles), as a function of the belt speed 
U, for H = 5.3 cm. The solid line is a fit to a linear variation 
of u) vs. U, consistent with a Hopf bifurcation. The solid 
circle marks the extrapolation to the critical value uic- 



based on symmetry and the assumption of a constant 
contact-point speed. 

We employ coordinates in the lab frame with x longi- 
tudinal to the belt, y transverse to it and z up. Ribe et 
al. [16] observed that the linear perturbation equations 
about a steady stretched catenary decouple into disjoint 
systems for perturbations in the plane (x, z) of the cate- 
nary and those out of the (x, z)-plane. Meandering is 
an out-of-plane perturbation at linear order. We observe 
that the reflectional symmetry of the catenary further 
implies that the nonlinear perturbation equations are in- 
variant to a change in the sign of all the out-of-plane 
variables. Hence if we perform a series expansion for the 
weakly nonlinear form of meandering near onset (^ ^ 1) 
then all the in-plane variables are even functions of A and 
all the out-of-plane variables are odd functions of A. In 
particular, the contact point {x{t),y(t)) between thread 
and belt must have the form 



x{t) 



Asm{ujt) + 0{A^), 

Xa + Bcos{2ujt + (j)) + 0{A^), (5) 

where B = 0{A^), uj = ujc+0{A'^), (p is a constant phase, 
and xo is the unperturbed displacement plus an 0{A^) 
correction. Thus, the absolute speed V{t) of the contact 
point is given by 



V' 



{x-Uf + v\ 
U'^ + 4BujU sm{2ujt + (f>) 
+ A^uj'^[l + cos{2ujt)]/2 + 0{A'^). (6) 



We now make the kinematic assumption that V is con- 
stant and always equal to the critical belt speed Uc- This 
can be motivated either by the empirical observation in 
Ref. P2] that the steady speed of coiling is remarkably 
close to the belt speed at marginal stability, or by an 



asymptotic argument that the amount of stretching in 
the bending boundary layer of a slender thread is negli- 
gible [18]. Setting V = Uc in Eqn.[6]gives 



U' 



U^ + A^u;y2 + 0iA^), 



(7) 



B = A^uj/{8U) + 0{A^), and q)=- + 0{A^). (8) 

We rearrange Eqn. [7] 

A^ = 2(C/2 - C/2)/c.2 « im!/u;l)[{Uc - U)m, (9) 
and compare with Eqn. |4]to predict 

ti^{cOc/2Ucf. (10) 

It is evident from Fig. [8] that this purely kinematic 
prediction of fi works remarkably well, particularly at 
larger heights. In effect, the buckling part of the thread 
behaves like an inextensible rope being played out onto 
the moving belt, though there is still extension in the 
falling part closer to the nozzle. The greater discrepancy 
between observation and kinematic prediction at lower 
heights is likely due to weaker scale separation between 
the falling and buckling parts of the thread. 

It would be interesting to examine the motion of an 
actual inextensible, flexible rope falling on the moving 
belt. Recently, the classic liquid rope-coiling instability 
has been compared to a real rope in just this way [19] . 
Presumably, the kinematic condition on the rope at the 
point of contact with the surface of the belt would have 
to be enforced by their frictional interaction. 

The above kinematic assumption also makes predic- 
tions for the amplitude B and phase ^ of a small, fre- 
quency 2uj, longitudinal oscillation. The existence of this 
oscillation can be discerned in the spacetime plot of me- 
andering viewed from the side, shown in Fig. Ma), which 
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FIG. 7: (Color online) Measurements of the Hopf frequency 
ujc (solid circles) as a function of the nozzle height H. The 
solid line is the theoretical result from the linear stability 
calculation of Ref. [16]. The reversals in the theoretical curve 
mark the multifrequency regime. 



is discussed in more detail below. The experimental reso- 
lution of B is however insufficient to make a quantitative 
comparison to theory. 

The form of Eqn. [5] is justified from consideration of 
the full amplitude equations, which are presented in Sec- 
tion Iml 



Beyond meandering: the figure eight and 
translated coiling states 



D. 



For H < 5.5 cm, the meandering state gives way at a 
non-hysteretic bifurcation to the figure eight state shown 
in Fig. |2];. The figure eight state is, in its turn, unsta- 
ble to the translated coiling state at lower H. This state 
links directly to pure fluid rope coiling a,t H — 0. It is 
difficult to extract quantitative information about these 
states from photos of the "stitch pattern" on the belt. 
However, by viewing the thread from the side, the x{t) 
and y{t) components of motion can be measured. We 
accomplished this by positioning a small mirror at 45° to 
the belt, so that a camera could image both the longitu- 
dinal and the transverse components at a point close to 
the location of maximum curvature of the thread, which 
was a few mm above the surface of the belt . Taking a sin- 
gle horizontal video line from such an image and plotting 
it in time gives spacetime images like those in Fig. [9] 

In these images, it is apparent that the meandering 
state consists of essentially a single frequency uj in y, with 
a small component at 2uj in x with amplitude B <^ A, as 
assumed in Eqn. |5] above. The bifurcation to the figure 
eight state clearly involves the excitation of a large 2uj 
mode in y, which mixes with the lo mode in x to produce 
the figure eight pattern. At the bifurcation to translated 



FIG. 8: (Color online) The nonlinear saturation coefficient fj., 
as a function of the nozzle height H. The solid circles are 
the values deduced from fits to Eqn. |4] Diamond symbols 
show the kinematic estimate fj, ~ {cJc/2Uc)'^ for each height, 
discussed in the text. The solid line is a cubic polynomial 
model of the kinematic estimates. 



coiling at lower H, the 2a; mode amplitude becomes small 
again, and the state is due to x and y oscillations at u) 
which are nearly 90° out of phase. 

These qualitative observations suggest that the mean- 
dering - figure eight - translated coiling series of transi- 
tions, and by extension all of the states of motion of the 
thread, could be understood in the framework of complex 
amplitude equations for modes with a few frequencies and 
their harmonics. These frequencies are probably closely 
related to the frequencies of pure coiling, including those 
in the multifrequency regime. 



III. AMPLITUDE EQUATIONS AND 
SYMMETRY 



We now discuss a more abstract approach to under- 
standing the "stitch patterns" through the machinery of 
equivariant bifurcation theory [20]. This theory offers a 
formal way to derive appropriate amplitude equations, 
starting from symmetry considerations alone. Such a 
qualitative approach offers a complementary perspective 
on the dynamics to that gained from fluid mechanical 
treatments. It also provides insight into the overall or- 
ganisation of the bifurcations between these highly sym- 
metric states. Given the rich phenomenology of this sys- 
tem, it seems wise to follow this route rather than at- 
tempt to account for all the states de novo from fluid 
mechanics alone. Although a rather large number of 
amplitudes may be required, this approach is still much 
simpler than the alternative. In principle, the full fluid 
mechanical description could be systematically projected 



onto a few modes, and the coefficients of the amplitude 
equations calculated. Alternatively, the coefficients could 
be deduced directly from experimental data of the kind 
shown in Fig. [9] 

We conjecture that the dynamics of the thread is close 
to the axisymmetric case of pure rope coiling, so that 
the motion of the belt can be treated as a perturbation. 
The moving belt breaks the rotational symmetry while 
preserving a reflection symmetry in a vertical plane con- 
taining the direction of travel of the belt. In bifurcation- 
theoretic terms, the onset of meandering from the steady 
catenary state is a Hopf bifurcation with weakly-broken 
0{2) symmetry [21]. We denote the frequency at the 
Hopf bifurcation by oj. To describe the secondary bifur- 
cation in which the figure eight states appear, a minimal 
description must include amplitude equations for a mode 
with frequency 2a;. We argue in the Appendix that the 
physical symmetries provide enough constraints to de- 
termine the form of the relevant amplitude equations. In 
section III A we discuss the effect of a further symmetry 
property which, while certainly not an exact symmetry 
of the physical problem, appears to help explain certain 
intriguing aspects of the dynamics of the viscous thread. 

We propose an abstracted description of the fluid me- 
chanical problem in which we focus on the position of 
the contact point of the viscous thread on the belt. We 
express the contact position of the viscous thread in the 
lab frame {x{t),y{t)) as the complex position coordinate 



Pit) 



x{t) + iy{t) 
+B+{t)e 



2iLut 



B-{ty 



-2iLjt 



(11) 



where A^, A_, i3_|_ and i?_ are the complex amplitudes 
of four rotating wave modes. By symmetry, the generic 
amplitude equations describing the dynamics, truncated 
at cubic order, take the form 



. = A+ (ei + iuji + ai\A+\^ + a2\A_\^ 




+bi\B+\^ + b2\B_\^) + siA_B+B_ 




+eb{eiA^ +PiA+B+ 




+qiA+B^ + p:iA^B+ + q^iA^B^), 


(12) 


= B+{e2 + iuj2 + ci\B+\^ + C2\B^\'^ 




+di\A+\'' + d2|A_|') + S2B^A+A^ 




+ebie2B^ +P2Ai + q2A\+piA+A^ 


-),(13) 



B, 



and two similar equations for A^ and i?_, given in the 
Appendix. The coefficients of the nonlinear terms in 
Eqns.[T2]and[T3]are in general complex- valued, and over- 
bars denote complex conjugates. The bifurcation param- 
eters ei and £2 are complicated unknown functions of the 
dimensionless groups identified above, or equivalently of 
the natural experimental control parameters U and H . 
ei{Uc,H) = at the onset of meandering, as described 
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FIG. 9: (Color online) Spacetime views of the side of the 
thread, with time running upward, a; is a coordinate lon- 
gitudinal to the belt, while y is transverse, (a) shows the 
meandering state, where almost all the motion is transverse. 

(b) shows the figure eight state in which two harmonics of the 
basic frequency are clearly visible in the longitudinal motion. 

(c) shows translated coiling, in which one main harmonic ap- 
pears, with a phase difference between x and y. 



above in section II B €2{Uc, H) remains negative and de- 
scribes the linear damping of the 2uj mode. The param- 
eter tb, a function of U with £{,([/ = 0) = 0, controls the 
strength of the symmetry-breaking effects of the belt on 
the thread. 



Eqns. [niand [13 
ical Eqns. [4[ andlS 



justify the use of the phcnomcnolog- 
as we discuss further below. A full 
analysis of the bifurcation structure of Eqns. |T2| and [T3| 
is mathematically feasible, but beyond the scope of this 
paper. In the remainder of this section, we first discuss 
the meandering instability, and then the transition from 
meandering to figure eights, in the context of these am- 
plitude equations. 

The steady catenary solution Aj^ — A^ ~ B+ — 
_B_ = exists for all parameter values and is stable if 
El ± efcRe(ei) < 0. Experimentally, there is a critical 
height Hq w 2.0 cm below which no rope coiling instabil- 
ity occurs for [/ = 0. It is therefore natural to speculate 
that, for small U and heights close to iJo, we could as- 
sume ei ex [H — Ho) and e^ oc U. These scalings in 
turn predict a straight-line marginal stability boundary 
for the onset of meandering near this endpoint. This 
regime is difficult to reach experimentally, however, due 
to the extremely slow time evolution and the thickening 
of the thread relative to its meandering amplitude. 

In the axisymmetric case U — 0, when 6^ = 0, there is 
an 0(2)-symmetric Hopf bifurcation a,t H = Hq in which 
two oscillatory branches appear |22| : a rotating wave cor- 
responding to coiling, and a standing wave correspond- 
ing to meandering. The coiling state persists under weak 
symmetry-breaking [21 and evolves smoothly into the 
translated coiling state A^ 7^ 0, A_ = B^ ~ i?_ = 0, 
observed for small U. In the axisymmetric case, the bi- 
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furcating standing wave has no preferred horizontal ori- 
entation. Under weak symmetry-breaking, two stand- 
ing waves out of this family of states survive. One of 
these states corresponds to oscillations where the thread 
lies entirely within the symmetry-plane of the belt at all 
times (a "longitudinal oscillation"), the other ("mean- 
dering") has the combination spatiotemporal symmetry 
K o Ttt (defined precisely in the Appendix) of reflection 
across the belt followed by waiting for a half-period of 
the oscillation. 

The figure eight state also has the symmetry kot-j^ and 
so, to simplify the further discussion of meandering and 
figure eights, we restrict our attention only to solutions 
that have the k o Tt^ symmetry, that is, which have the 
form A+ = -A_ = A and i?+ = B_ = B. We refer to 
the invariant subspace {A, —A, B, B) as S; importantly, 
S remains invariant even if higher-order terms are added 
to the amplitude equations. 

Within S, the amplitude equations simplify to become 

A = {ei+iuJi-ei,ei)A + iei,aiAB + l3iA\A\'^ + I32A\B\'^, 

(14) 

B ^ {62 + iuj2 + €1,62) B + ieba2A'^ + ^3B\B\^ + PiB\A\'^. 

(15) 

The coefficients ei, 62, ai,a2,/?i, • ■ • ,/34 in Eqns. [14] 
and [15] c an be deduced by direct comparison with 
Eqns. |12| and [13) Near the initial meandering instabil- 
ity ei — ef,Re(ei) = 0, and we may eliminate B adia- 
batically, or via a centre manifold reduction, so that to 
leading order B = —i€ba2A^/{e2 + iuj2 + £662) -t- 0{A^) 
and hence 



to ensure the existence of an instability of meandering 
to the 2uj mode that produces the figure eight pattern. 
At €1 = Re(/3i)e2/Re(/34) there is a pitchfork bifurcation 
at which meandering loses stability to "mixed-modes" 
involving the 2lu mode. For eb ^ the pitchfork bifur- 
cation becomes imperfect, with a continuous, but rapid 
and non-hysteretic, increase in the amplitude of the 2uj 
mode in the vicinity of ei = Re(/3i)e2/Re(/34). 

Thus, the sequence of transitions observed experimen- 
tally has a natural interpretation in terms of the mode 
interaction described above. 

Moreover, the experimental evidence indicates a dis- 
continuous but not hysteretic transition from meander- 
ing to figure eight; from the theory one would expect 
either that the transition is indeed abrupt (for example a 
subcritical bifurcation) and accompanied by measurable 
hysteresis, or that the transition is non-hysteretic but 
continuous. The resolution of the exact nature of the 
transition from meandering to the figure eight state will 
require better controlled experiments to sort out. Direct 
measurement of the amplitudes and phases of the lo and 
2uj components of the pattern should enable this intrigu- 
ing discrepancy to be resolved. 



(ei -f- iuji — ebei)A + 



ef^aia2 



£2 + iUJ2 + £662 



Pi]A\A\' 



(up to 0{A^)), which is the usual amplitude equation for 
a Hopf bifurcation, justifying the use of Eqn. [4J 

Writing A = Eie^^i and B = i?2e*^2 we find that, due 
to the time-translation symmetry, Eqns. [141 and |15| reduce 
further, to a three-dimensional set of ODEs for Ri, R2 
and the relative phase X = ^2 — 20i . There are regions of 
parameter space in which time-periodic solutions with a 
constant relative phase x do not exist. In other regions 
of parameter space, there will typically be two possible 
constant values x = X±- ^^ general, within S we recover 
the form 



y{t) ~ 2Rism.ujt, 

x{t) = 2i?2Cos(2wi + x), 



(16) 
(17) 



for the location of the tip of the thread, which confirms 
Eqn. [5] 

If ei > and €2 < (as expected for the onset of 
meandering when H < 6.0 cm), and if e^ — 0, we 
would expect stable meandering to exist for < ei < 
Re(/3i)e2/Re(/34), choosing Re(/3i) < to ensure that 
the onset of meandering is supercritical, and Re(/34) > 



A. An additional model symmetry 



The amplitude equations, discussed above, do not pre- 
cisely determine the relative phase x- However, experi- 
mental results (Fig. [2[as well as those of Ref. 15j) indi- 
cate that the values x — '^ 1"^ ( the "bunched-up mean- 
dering" state discussed in Ref. [TS]) and near x = — 7r/2 
(figure eight) are consistently selected. It is natural to 
ask whether there is a reason for this, either of an ab- 
stract or a fiuid-mechanical nature. 

As far as the abstract model problem goes, it turns out 
that there is indeed an additional symmetry that implies 
X = ±7r/2. This symmetry (labelled h in the Appendix) 
is a 180° rotation of the (x, y) plane followed by reversal 
of the direction of time. Clearly this is not a symmetry of 
the original physical problem: in no sense are we suggest- 
ing that falling fluid in the thread is related to a viscous 
thread of rising fluid! But in the bifurcation-theoretic 
model, we ignore the vertical direction, and the time- 
reversal applied to the equilibrium state results only in 
the reversal of the direction of travel of the belt. Even so, 
we would expect h to be only an approximate symmetry 
since the initial steady catenary state is not symmetric 
under a 180° rotation. The consequences of the extra 
symmetry h would be to force the imaginary parts of all 
coefficients in Eqns. [14] and [15] to be real. In this case, 
the ODEs for R\ , R2 and the relative phase x = ^2 — 26'i 
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become 
R2 



(ei - ebei)Ri - ebaiRiR2smx + I3iRI + /32i?ii?|, 

(18) 

(£2 + e6e2)i?2 + eba2Rl sinx + /33^2 + P4R2RI, 

(19) 



X = 



£(, I — 2aii?2 I cosx- 



(20) 



Hence, from Eqn. [20] time-periodic solutions with a con- 
stant relative phase exist only for x = ±7r/2. If the coeffi- 
cients in Eqns. 14 and 15 have (small) nonzero imaginary 



parts then the two time-periodic states persist but now 
only have % « ±7r/2. There is some experimental evi- 
dence that the bunched-up meandering state observed in 
Ref. |15| has x ~ '^/'^i without having exact equality. 



More generally, Eqns. 18 



20 are equivalent to the am- 
plitude equations of the well-studied 1 : 2 steady-state 
mode interaction |521 IMl US [55] . The meandering and 
figure eight states correspond to "mixed-mode" equilib- 
ria in this language, since both states contain non-zero 
amounts of both the u; and 2io modes. However, the 
analysis in Refs. [23J [231 HSl I2S] concentrates on cases 
where the quadratic terms have order unity coefficients 
and strongly influence the dynamics near Ri — R2 — 0. 
In contrast, in the present case, we are interested in the 
regime e^ <C 1 where we see from Eqns. [18] - [20] that the 
quadratic terms flx the relative phase x but only weakly 
affect the amplitudes Ri and i?2- 

Having derived the form of the amplitude equations 
on symmetry grounds, an obvious next step would be 
to compute or measure the undetermined coefficients in 
them. It is straightforward to extend Eqns . \T2\ and \T3\ to 
account for higher frequencies, although the number of 
terms, even at cubic order, grows quite rapidly. Project- 
ing the full dynamics of the thread onto these equations 
serves to simplify the problem to finding the coefficients 
and how they depend on the nondimensional experimen- 
tal parameters H{g/v'^Y/^, U{gH^/v) and Hi, H2 and 
Ha given by Eqns. [l][3] In principle, the entire state dia- 
gram shown in Fig. Hi may be understood in this way. 



IV. CONCLUSIONS 

In this paper we have reported a detailed experimental 
survey of the "stitch patterns" observed in a more precise 
version of the "fluid-mechanical sewing machine" |15| . 
We studied, in particular, the onset of the meandering 
state and measured its critical belt speed, Hopf frequency 
and amplitude saturation. Our theoretical understand- 
ing and interpretation of this instability of the catenary 
thread comes from three sources: a detailed numerical 
solution of the linear stability problem |16j . an ad hoc 
argument which used the assumption that the contact 
point of the thread with the belt has a constant velocity, 
and the construction of the full amplitude equations on 
symmetry grounds. 



We found that the onset of meandering is very well 
described as a Hopf bifurcation with the critical param- 
eters Uc and tOc accurately predicted by a linear stabil- 
ity analysis of the fluid equations [16J . Surprisingly, the 
nonlinear saturation of the meandering motion could be 
mostly accounted for by simply assuming that the fluid 
thread behaves as an inextensible rope being played out 
onto the belt. This is the case despite the considerable 
buckling and stretching of the real thread. Finally, the 
full amplitude equations, derived from symmetry con- 
siderations independent of the fluid mechanical details, 
provide a systematic way of organizing and interpreting 
the otherwise rather complex motion of the thread. 

From the full amplitude equations we were able to 
explain, at least in outline, the existence of the figure 
eight bubble that occurs over a substantial part of the 
state diagram for smaller H, as seen in Fig. [3] It seems 
likely that this state is a manifestation of a 1 : 2 steady- 
state mode interaction [23] IMJUSl US] ■ In general, such 
a mode interaction is well known to generate complex 
dynamics in certain cases, such as when ai and a2 in 
Eqns. [18] - [20] have opposite signs. In this case, one ex- 
pects quasiperiodic oscillations and robust heteroclinic 
cycling. Physically, such dynamical complexity would be 
most easily observed in the intermittent evolution of the 
relative phase variable x on a timescale that is slow com- 
pared to the basic meandering oscillation period. Future, 
improved experiments could look for such subtle effects. 

An interesting question that needs to be addressed 
more fully, both in experiment and theory, is the vari- 
ation of the oscillation frequency as the belt speed is 
reduced. Ribe ei. al. |16j compared the frequency at 
the onset of meandering to the frequency of finite am- 
plitude "pure" coiling (when C/ = 0) at the same height 
H and found excellent agreement. Our results, however, 
indicate that the meandering frequency decreases sub- 
stantially from that at onset as the belt speed decreases, 
as shown in Fig. [6] This decrease is exactly what is ex- 
pected from the linear amplitude dependence of the fre- 
quency following a Hopf bifurcation. But how can this be 
reconciled with the requirement that the frequency must 
also agree with the pure coiling frequency again when U 
is ultimately reduced to zero? Is it the case that the fre- 
quency reaches a local minimum as U is decreased, before 
increasing continuously back towards pure coiling, or are 
there discontinuities? What is the effect of the interven- 
ing figure eight state on the fundamental frequency? 

For larger values of H, the amplitude equation model 
would have to be extended to allow for modes with the 
additional frequencies that arise from the new pendulum 
motions of the thread. Using data gathered from side- 
view images like that shown in Fig. [9] amplitudes and 
phases of the constituent motions of the thread could be 
acquired and interpreted within the framework of the pre- 
dictions of the amplitude equations. Future work using 
this approach will require a more complete investigation 
of the bifurcation structure of multifrequency amplitude 
equations with weakly broken 0(2) symmetry. 
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APPENDIX: RESONANT HOPE BIFURCATION 
WITH WEAKLY BROKEN 0(2) SYMMETRY 

In this Appendix, we briefly sketch the derivation of 
the abstract amplitude equations, Eqns. 12 and [T3| de- 



scribing the dynamics for heights H towards the lower 
end of the range investigated experimentally. We em- 
phasize that this discussion proceeds in a manner inde- 
pendent of the fluid-mechanical details of the flow, and, 
although it provides some support for our experimental 
results and simple explanatory ansatz, it remains a mod- 
elling approach. 

The system is taken to be a weak perturbation of the 
axisymmetric case that corresponds to U = 0. We first 
compute the amplitude equations for the axisymmetric 
case and then discuss how the weak symmetry-breaking 
provided by the belt admits additional coupling terms. 

In the axisymmetric case, the symmetry group con- 
sists of planar rotations about the origin, and reflections 
in vertical planes containing the origin: this is the or- 
thogonal group 0(2). The action of 0{2) on the plane 
E^ containing the belt is generated by the elements k (a 
reflection in the plane containing the belt) and pe (ro- 
tation through angle 9). Identifying {x,y) E M^ with 
p = X + iy G C we have 



k{p) = p, 



Pe{p) 



e^'p. 



In addition, the system has no preferred origin in time, 
and hence the time-periodic solutions we investigate have 
no preferred phase. This additional symmetry r^ : i — > 
t + (p/uj generates a group S^ corresponding to phase 
shifts of the ampHtudes. We suppose that the coihng 
motion can be described as a sum of four oscillatory am- 
plitudes for a complex position coordinate, as in Eqn. |ll[ 
Let A = {A^,A^,B+,B-) denote the vector of ampli- 
tudes. The action of the full symmetry group 0{2) x S^ 
on A is inherited from the action on M^ of the generators 
{K,pg), and the phase-shift symmetry r^: 



that are symmetric under Eqns. [X?T] - |A.3| obey the equiv- 
ariance condition 7F(A) = F(7A), for all group elements 
7 £ 0{2) X S^. The equivariance condition is a mathe- 
matical statement of the intuitive notion that whenever 
A(t) is a solution of the amplitude equations, we require 
7A(t) also to be a solution trajectory. We can now deter- 
mine the form of the most general amplitude equations by 
computing invariant and equivariant polynomials. Gen- 
eral methods for this exist [^U] but are too cumbersome 
to apply here, particularly since we intend only to com- 
pute terms at the first few orders. As usual in problems 
of this type, the terms \A+\^, |A_|2, jB+p and |B_|2 
(and sums and products of these) are invariant. To find 
other invariants we adopt the "brute force" approach of 
looking for an invariant term A'^A'^B'^B'^i with integer 
powers, using the convention that a negative power de- 
notes a positive power of the complex conjugate variable. 
Invariance under pg and r^ implies a + b + c + d = and 
a — b + 2{c — d) = 0. We define the order of the in- 
variant term to be \a\ + \b\ + \c\ + \d\. The invariants 
of order 2 are \A+\^, |A_|2, \B+\^ and |B_|2. It is eas- 
ily seen that the only invariant of order 4 that is not a 
product of order 2 invariants is A+A-B+B-. It follows 
that the terms up to cubic or der i n the first component 
ii =Fi{A+,A^,B+,B^) of ([aHI will be 



k{A+,A^,B+,B^) 
pe{A+,A^,B+,B^) 
T^[A+,A^,B+,B^) = {c'^A 



The required amplitude equations 
A = F(A), 



,e2^"^B+,e-2*'^B_). 
(A.3) 



(A.4) 



A+\B, 



A+\A, 
AAB^ 



A+\A.\^ 
A_B+B_ 



By now applying the reflection symmetry k, 
we may deduce the second amplitude equation 



F2{A+,A_ 



,B+,B^ 



-). Note that n implies that the 



coefficients in F2 are the complex conjugates of those in 
Fi. 

Similarly we deduce the order 2 and order 4 invariants 
for F3 and then use k to compute F4 from F3. The result- 
ing amplitude equations, Eqns. 12 and 13 with e^, = 0, 



describe the axisymmetric interaction of the ui and 2uj 
modes. 

The introduction of a (slowly-moving) horizontal belt 
changes the symmetry of the problem in two ways. The 
first of these is clear: it breaks the rotation symmetry pg. 
To account for this effect we compute low-order terms 
that are not symmetric under pg but which do respect 
the remaining symmetry r^. The eight new invariants 
are: 



{A^,A+,B^,B+), (A.l) 

e^'>iA+,A^,B+,B^), (A.2) 



A+A^B+, AtB+, 
AlB+, A+A^B^, A^_B^, AlB^ 



A+A^, B+B 
\lB+ 



The new quadratic invariants give rise to the new linear 
terms in Eqns . [T2| and [T3] and so split the symmetric Hopf 
bifurcation at ei = into two generic Hopf bifurcations 
at ei ± ef,Re(ei) = 0. The new cubic invariants generate 
the new quadratic terms. The amplitude equations for 
the four amplitudes are then given by Eqn. 12 and 13 
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along with 



i- = A_ (ei -iwi + ai|A_p + a2|A+p 

+ebieiA++piA_B_ 
+qiA^B++p3B^A+ + q3A+B+), (A.5) 



B. 



= S_(e2-ZCJ2 + Ci|B_|2+C2|B+|2 

+Jl|A_|2 + J2IA+P) + S25+^+A_ 

+eb(e2-B+ +P2^+ + 92^- +P4^+A_)(A.6) 



Some features of the phase relationship between the lu 
and 2a; modes suggest the introduction of an additional 
symmetry, denoted h: rotating the apparatus by 180° 
around the vertical z-axis and reversing the direction of 
travel of the belt (i.e. reversing the direction of time) 
leaves the basic state p = unchanged. For the bifurca- 
tion problem constructed in this Appendix, h is an inde- 
pendent symmetry that introduces a new constraint. But 
h is clearly not an exact symmetry of the original fluid- 
mechanical problem where the fluid in the viscous thread 
falls under gravity, and the initial equilibrium state is a 
catenary rather than an exactly vertical thread. So, at 
best, the effect of h is felt in some approximate way by 
amplitude equations describing the original physical sys- 
tem. 



The new symmetry h{x,y,t) 
mode amplitudes as 



-{x,y,t) acts on the 



h{A+,A_,B+.B_) = -{A_,A+,B_,B+). (A.7) 



So the composite operation h o k sends each amplitude 
to the negative of its complex conjugate; equi variance 
of Eqn. |A.4| under h o k implies that odd order terms 
would have real coefficients and even order terms would 
have purely imaginary coefficients. More realistically, we 
might hope that the coefficients have non-zero but small 
(compared to their modulus) imaginary or real parts at 
odd and even orders respectively. 



The analysis of Eqns. [nj [13] |A3] and |A]6] is simplified 
by the existence of invariant subspaces for the dynamics - 
technically, the fixed-point subspaces of group elements. 
The most relevant of these for this paper consists of those 
points left unchanged by the combined operation k o r^r 
(reflection in the plane of the belt followed by a time- 
shift of half an oscillation period) : it is therefore useful 
to define the invariant subspace 

S = Fix{KOT^) = {A+,~A+,B+,B+), 



which contains both the meandering and figure eight 
states. 
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